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Clark  (1976)  presents^-a  simple  method  by  wliich  to  calcu- 
late the  change  in  sea  level  caused  by  tlie  gravitational  attrac- 
tion of  an  ice  cap.  In  the  same  paper  he  etlso  presents  some 
results  of  a more  complex  numerical  solution  of  the  problem, 
and  the  two  calculations  are  in  disagreement.  ^ In  fact;  his 
more  powerful  solution  is  correct.  His  simple  calculation  con- 

A 

tains  an  error.  It  may  be  useful  to  pointT”out  the  reason  for 
this,  lest  doubt  be  cast  on  his  numerical  method  and  the  many 
interesting  consequences  of  his  thesis. 

When  a point  mass  M is  created  on  the  surface  of  a per- 
fectly ri^gid  earth,  Clark  concludes  that  the  change  in  sea 
level  e degrees  away  from  the  mass  is 

S(e)  = GM  csc/Qv 
2ag  ' 2 > 

where  G is  the  gravitational  constant,  a is  the  radius  of  the 
earth,  and  g is  the  acceleration  of  gravity  at  the  earth’s  sur- 
face. This  might  better  be  written  as 

S (0)  = ^ csc/9  V 
2^2' 

where  x is  the  ratio  of  M to  the  mass  of  the  earth.  This  equa- 
tion is  obviously  incorrect,  as  s is  always  positive,  leading 
to  a rise  in  sea  level  everywhere.  This  cannot  happen  if  the 


iTiasu  of  the  ocean  is  conserved,  as  Clcirk  demands.  The  error 
enters  Cltirk's  derivation,  described  in  his  Appendix  1,  in  the 
equation  (s)  ~To'  that  is,  tliat  the  gravitational  potential 
on  the  sea  surface  cquipotential  after  the  creation  of  the 
mass  is  equal  to  that  on  the  sea  surface  equipotential  before 
the  creations  of  the  mass.  If  the  mass  of  the  ocean  is  to  be 


held  constant  this  will  not  be  true.  The  two  sea  surfaces  will, 
of  course,  be  equipotential  surfaces,  biit  their  potential  values 
(relative  to  a zero  valiie  at  inf inite  radius)  v;ill  in  general 
differ. 

Clark's  equation  does  give  the  correct  shape  for  an  equi- 
potential surface  near  the  sea  level  surface,  and  one  can  use 
it  to  reach  the  correct  solution  by  making  tv;o  assumptions: 

1)  In  the  region  of  sea  level,  the  gravitational  gradi- 
ent is  of  equal  magnitude  every\i7here. 

2)  In  terms  of  the  varia}>le  0,  oceans  and  continents 
are  uniformly  distributed  over  the  earth. 
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where  C is  a constant  whose  value  can  be  determined  by  requir- 


ing no  not  change  in  the  v'Olume  ol  the  ocean » 
rTT 
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Tiiis  leads  to 


sin  e 

d-e  = 

0. 

= ax 

f CSC 

(£)-  2 > 

- 

2 

1 

'2'  j 

Thus,  gravitational  forces  cause  sea  level  to  rise  over  the  area 
within  60°  of  the  mass  and  to  fall  over  the  rest  of  the  sphere. 
One  can  get  an  idea  of  the  significance  of  the  second  assump- 
tion by  noting  that  over  a wide  range  of  ocean  distributions 
(from  the  near  hemisphere  all  ocean,  and  the  far  hemisphere  ail 

continent  to  the  opposite  case),  C varies  from  -2.85  a^c  to 

^ 2 

-1.17 

2 ' 

Clark  uses  his  method  to  calculate  the  gravitation  compo- 
nent of  sea  level  change  near  Greenland  due  to  the  melting  of 
Greenland  ice  cap  since  9000  years  B.P.  He  finds  a 27  m compo- 
nent. If  one  uses  the  more  correct  equation,  one  finds  a 25  m 
component.  The  difference  between  the  two  results  is  small  be- 
cause the  mass  of  ice  involved  is  small  (ax  g 1.07  m)  . A muc); 

2 

larger  mass  of  ice,  however,  can  lead  to  serious  discrdpenci es . 
Using  a more  complex  numerical  metliod,  Clark  considers  the 
melting  of  the  Laurentide  and  Fennoscandian  ice  caps.  This  nu- 
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merical  npproacOi  includGa  the  offcctn  of  the  lateral  extent  of 
the  ice  caps  and  the  self-gravitation  of  the  oceans.  If  the 
inciting  raises  the  average  sea  level  by  85  m (ax  = 16.4iii), 

2 i 

Clark's  more  siitiplc  method  predicts  sea  level  rise  of  69  m 
in  the  South  Pacific  and  52  m about  60°  from  the  ice  cap  cen- 
ter. This  disagrees  with  the  results  of  his  numerical  analy- 
sis, as  shown  in  his  Figure  2.  The  method  derived  here,  how- 
ever, predicts  a rise  of  101  m in  the  South  Pacific  and  85  m 
about  60°  from  the  ice  cap  center,  in  complete  agreement  with 
Clark's  numerical  solution.  This  indicates  that  Clark's  numer- 
ical approach  is  correct,  in  spite  of  the  error  in  his  simple 
calculation.  If  Clark  has  not  included  the  land-sea  distribu- 
tion in  his  solution,  however,  and  he  does  not  state  that  he 
has,  he  should  do  so  before  continuing  with  his  stated  plan  of 
refining  his  model  to  look  at  small  fluctuations  in  ice  sheets. 

^•Thile  the  correction  v;ill  have  little  effect  on  the  size  of 
small  fluctuations  in  local  sea  level,  it  will  noticably  ef- 
fect large  changes.  including  the  correction  v;ill  be  generally 
useful  and  v;ill  add  to  the  value  of  the  result. 
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